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Tunneling of electrons traversing a few-electron quantum 
dot is strongly influenced by the Coulomb interaction lead- 
ing to Coulomb blockade effects and single-electron tunneling. 
We present calculations which demonstrate that correlations 
between the electrons cause a strong suppression of most of 
the energetically allowed tunneling processes involving excited 
dot states. The excitation of center-of-mass modes, in con- 
trast, is unaffected by the Coulomb interaction. Therefore, 
channels connected to these modes dominate the excitation 
spectra in transport measurements. 



Transport measurements are a highly sensitive tool for 
the investigation of the electronic structure of semicon- 
ductor quantum dots. While the differential conductance 
of such a system in the linear transport regime is domi- 
nated by the classical Coulomb blockade effect Q , quan- 
tum mechanics leaves its fingerprints in magnetic field de- 
pendent addition spectra ||^,^ and in excitation spectra 
which can be obtained from conductance measurements 
in the finite drain-source voltage regime |^-^ . 

Due to a small electronic effective mass, quantum dots 
in semiconductor heterostructures exhibit a discrete level 
spectrum, easily discernible at low temperatures. Since 
the charge of these quantum dots can be controlled at 
the single-electron level, they are often called artificial 
atoms. The lateral confinement potential which binds the 
electrons to the quantum dot is typically created by spa- 
tially extended charge distributions. Therefore, it obeys 
a parabolic dependence on the distance from the cen- 
ter of the system (oc r^) rather than the f /r-dependence 
characteristic of the core potential of natural atoms 0. 
In analogy to the optical spectroscopy on conventional 
atoms the first attempts to study the excitation spectrum 
of these artificial atoms have been made by far-infrared 
(FIR) spectroscopy However, the long wavelength 

of the far-infrared radiation together with the parabolic 
confinement potential prohibits most transitions into ex- 
cited states by strict dipole selection rules |fj,|l^ . 

Transport spectroscopy on the other hand is expected 
not to suffer from these restrictions, but to allow a unique 
full spectroscopy of the level structure in these artifi- 
cial atoms. Nevertheless, experiments on few-electron 
quantum dots exhibit only a sparse excitation spectrum 



with relatively large level spacings ||]. These measure- 
ments suggest that the dominant resonances are due to 
a constant level spacing, independent of the number of 
electrons in the quantum dot. This behavior, which 
could be understood for a non-interacting electron sys- 
tem, is indeed remarkable, since the Coulomb interaction 
strongly influences the energy spectrum of the quantum 
dot and leads to small level spacings and strong correla- 
tions |T|,|I|. 

Thus, the experimentally observed large level spacing 
together with the appearance of a characteristic excita- 
tion energy demands an explanation. We present calcu- 
lations which show that even in transport experiments 
the excitation spectra are dominated by center-of-mass 
modes, similar to the situation in FIR spectroscopy. This 
gives rise to the constant level spacing with the charac- 
teristic single-particle excitation energy. This situation 
is caused by a suppression of most other transitions due 
to strong correlations between the electrons. 

We study the transport properties of a two-dimensional 
quantum dot coupled to two reservoirs by tunneling bar- 
riers in the conventional tunneling Hamiltonian approach 
| p5[ . The Hamiltonian of the system is given by a sum 
of the dot Hamiltonian Ho, a reservoir Hamiltonian Hn, 
and a tunneling term Ht- The tunneling Hamiltonian de- 
scribes the transfer of electrons from the reservoir, where 
the Coulomb interaction between the electrons is effec- 
tively screened (metallic regime), to the quantum dot, 
where interactions are most important. Our results are 
based on features of the few-particle eigenstates of the 
quantum dot which we obtain from a numerical diago- 
nalization of the dot Hamilton operator 



Vr,: 



i^dXd^dn'dm' ■ (1) 



It describes interacting electrons in a parabolic confine- 
ment potential subjected to a perpendicular magnetic 
field B. The single-particle energies e„ are connected 
with left and right circularly polarized oscillator eigen- 
modes, 



hn+ {N+ + 1/2) + nn^ (n^ + 1/2) 



(2) 



with n± = (-y/^ri^+cjj ±ct'c)/2. r^o characterizes the 
confining potential and coc = eB/imc) is the cyclotron 
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frequency. A crucial feature of the parabolic confine- 
ment potential is the separation of the dot Hamiltonian 
into a relative and a center-of-mass (CM) part Jl6| , p7t , 
where the spectrum of the CM-Hamiltonian is identical 
to the single-particle spectrum (|^) and independent of 
the number of electrons in the dot. It is this separation 
which prohibits in FIR spectroscopy the observation of 
the spectral fine structure of these artificial atoms, and 
makes visible only the CM-excitations. The Coulomb in- 
teraction between electrons enters only the Hamiltonian 
of the internal degrees of freedom (relative part). We ob- 
tain the matrix elements for the interaction, Vnmn'm', in 
closed form for the usual l/{sr) Coulomb potential with 
the dielectric constant e. 

Regarding only sequential tunneling, the current 
through the quantum dot is given by |15|| 

I^-eJ2 Taa' [P{N, a) + P{N - I, a')] 

a.a' 

X [f{AEaa' - mO - ./(A^aa' - Mr)] (3) 

where the resonant energy AEaa' = E{N,a) — E{N — 
1, a') is the difference between the energy of a iV-particle 
state a and a {N ~ l)-particle state a' . The Fermi-Dirac 
distribution function / characterizes the occupation of 
electron levels in the left (electro-chemical potential /i/) 
and right (/j.^) reservoir. The probability P{N, a) to find 
the quantum dot in the TV-particle state a will deviate 
from its equilibrium value for a given drain-source volt- 
age (/i; — /Zr)/e. Its dependence on the tunneling rate 
Taa' is well described by kinetic equations ||l5| , [l9| and 
leads to the blocking of conducting channels and nega- 
tive differential conductance p|, |l5| , pO| . 

In a "random phase" approximation, i.e. neglecting 
phase correlations between the initial state of the elec- 
tron in the lead and its final state in the quantum dot, 
the tunneling rate Taa' factorizes into an effective tun- 
neling rate for a non-interacting electron traversing the 
barrier, 7, and an overlap or spectral weight matrix ele- 
ment. ||T^ (For simplicity, we assume a constant effective 
tunneling rate and neglect the dependence of 7 on the 
single-electron dot states.) The overlap matrix element 
is given by 

Y,\{N,a\dt\{N-l),a')\\ (4) 

n 

This quantity describes to what extent the compound 
state, built by an incoming electron and the {N — 1)- 
electron state a' in the dot, overlaps with the iV-electron 
dot state a. While for an uncorrelated electron system 
this overlap evaluates to unity, correlations reduce it con- 
siderably [TslJl^ ]. Due to the summation over all single- 
particle dot states with equal weight in cq. (^ this quan- 
tity is insensitive to any features of the incoming electron. 
It only reflects the correlations in the states a and a' . In 
the following we will concentrate on these overlap matrix 



elements, as they predominantly determine the transport 
through the quantum dot (for levels separated by more 
than /cT, as discussed for example in Ref. JTst . In exper- 
iments, kT « 50mK « AfieV). Our numerical results are 
obtained for quantum dots with up to A'^ = 3 electrons. 

Fig. |l] shows the spectra of a quantum dot occupied by 
one, two or three electrons (QD-hydrogen, QD-hclium, or 
QD-lithium) at zero magnetic field. Due to the Coulomb 
interaction the level spacing of the few-particle spectra 
is much smaller than the single-electron level spacing, 
which is given by the confining energy hU,Q = 2 meV (this 
value is characteristic of experimental geometries, such 
as the system in Ref. Q). Enforced by energy conserva- 
tion, tunneling through the quantum dot is only allowed 
when the energy of the incoming particle matches the en- 
ergy difference between a QD-helium (lithium) eigenstate 
and a QD-hydrogen (helium) eigenstate. The occurrence 
of all energetically allowed transitions would clearly re- 
sult in a dense excitation spectrum. That contradicts 
recent experiments [|. Thus, it is not only the energy 
requirement which determines whether an electron tun- 
nels through the quantum dot. A mechanism to select 
certain transitions is provided by the overlap matrix ele- 
ments, eq. (||). 

One of these selection rules is evident from spin conser- 
vation pO| . For our few-electron quantum dots they for- 
bid transitions between a spin singlet state of QD-helium 
and a spin quartet state of QD-lithium. This has al- 
ready consequences for the usual Coulomb blockade oscil- 
lations, as transitions between the ground states of QD- 
hclium and QD-lithium (G-G-transition) in a perpendic- 
ularly applied magnetic field are considered. The total 
spin of the ground state of these few-electron systems os- 
cillates between its maximum (triplet, respectively quar- 
tet) and its minimum value (singlet, respectively doublet) 
when the magnetic field is increased Fig. |^ shows 

the ground state energy of QD-helium and QD-lithium 
as a function of the magnetic field B. The total spin is 
indicated for the different field ranges. For most values 
of the magnetic field G-G-transitions are possible. How- 
ever, as indicated in the figure, two field ranges occur 
where the ground state of QD-helium is a singlet while 
that of QD-lithium is a quartet: the G-G-transition is 
blocked. Thus, a drastic reduction of the conductance in 
the linear response regime should be observed when the 
magnetic field is tuned into these ranges. 

However, spin selection rules alone are not sufficient to 
explain the sparseness of the observed excitation spectra. 
Apart from an angular momentum conservation in the 
dot, the effect of which is annihilated by connecting the 
dot laterally to the leads, the overlap matrix elements, 
eq. (^), provide no further strict selection rules. Never- 
theless, very small spectral weights strongly reduce the 
probability for the occurrence of the corresponding tran- 
sitions. This gives rise to quasi- selection rules. Fig. ^ 
shows the values of the overlap matrix elements for tran- 
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sitions from the QD-helium ground state to all states of 
QD-lithium at -B = 0. They correspond to the most im- 
portant processes, if the quantum dot is in its ground 
state before the subsequent tunneling process. Since the 
2-particle ground state is a singlet state, only transitions 
into doublet states are possible. Transitions into degen- 
erate states are subsumed under the same peak, since all 
corresponding overlap elements contribute accordingly to 
the tunneling rate at a given energy AE. This applies 
especially to the double degeneracy with respect to the 
z-component of the doublet spin, giving rise to a factor of 
two in the peak heights. At low energies, below 15 meV, 
the figure exhibits a few regularly-spaced dominant peaks 
with several relatively smaller ones between them. This 
energy regime coincides with the experimental situation 
where only one electron tunnels at a time |^]. 

The first peak at 9.7 meV corresponds to the G-G-tran- 
sition. The B = states are, apart from the spin degen- 
eracy, degenerate with respect to the sign of the angular 
momentum. Therefore, the 3-particle ground state (with 
an orbital angular momentum of h) is four-fold degener- 
ate. Each of these states contributes equally to the peak 
shown in Fig. p[ Thus, the individual overlap matrix el- 
ements (eq. (0)) are much smaller than unity, indicating 
the strong correlations in the few-particle states. Nev- 
ertheless, compared with the other spectral weights the 
individual overlap matrix elements of the G-G-transition 
are the largest. Therefore, these transitions occur with 
the highest probability. They are followed by a transition 
which is subsumed in the peak appearing at 11.7 meV, 
i.e., at an excitation energy of 2 meV = hQa above the 
G-G-transition. This transition corresponds to a state of 
the 3-particle system with one quantum of CM-excitation 
in addition to the ground state energy. As a conse- 
quence of the separation of the CM- and relative mo- 
tion, each transition within the internal degrees of free- 
dom is accompanied by independent excitations of CM- 
modes. Since the excitation spectrum of the CM-motion 
is equidistant, these transitions give rise to equidistant 
conductance peak replicas (see arrows in Fig. I). More- 
over, at zero magnetic field the CM-energy spectrum is 
highly degenerate. Therefore, the overlap matrix ele- 
ments of several transitions contribute to the tunneling 
probability for a given resonant energy. It is this preva- 
lence of the CM-excitations which provides an explana- 
tion for the experimental observation of a characteristic 
level spacing independent of the number of particles in 
the quantum dot: The CM-spectrum is independent of 
the number of particles in a parabolic quantum dot and 
equals the single-electron spectrum ||l2| . [Notice that the 
small thermal broadening of the Fermi level in experi- 
ments, typically 100 mK, will not blur the different CM- 
excitations, occurring with much larger spacings, 2 meV 
here.] 

It is clear that in addition to the G-G-transition, tun- 
neling via other internal excitation states is possible, with 



their corresponding CM-excitation replicas (see second 
pair of replicas in Fig. ^). In some cases (e.g. in the 
transition from the single-particle ground state to the 
first-excited 2-particle state, not shown here) the over- 
lap matrix elements can even be larger than those of the 
G-G-transition. Due to their large degeneracy, transi- 
tions which end in a maximum-spin state produce very 
pronounced series of replicas, much stronger than the 
ones shown in Fig. ^. The occurrence of more than one 
group of dominant replicas gives rise to smaller level spac- 
ings in the spectra. Both effects have been observed in 
experiment for different iV-states in the dot. However, 
we found in all cases that only few groups of dominant 
replica occurred, while most excitations of internal de- 
grees of freedom are suppressed in tunneling. In our 
model, with a constant tunneling rate 7, this suppres- 
sion is solely due to correlations of the few-particle wave 
functions. These strong built-in correlation effects lead 
to nearly- forbidden transitions. Notice that the effective 
suppression of internally-excited states becomes less pro- 
nounced in the high energy range, as indicated by a more 
even distribution of high-overlap peaks in Fig. ^. 

Although our calculations only include systems with 
up to three electrons, general conclusions can be drawn 
for quantum dots occupied by more electrons. As long as 
correlations strongly influence the low energy excitations, 
a strong reduction of most overlap matrix elements oc- 
curs, leading to quasi-selection rules. These will suppress 
most transitions involving excitations of the internal de- 
grees of freedom in the rather dense spectrum. Since 
the CM-motion remains unaffected, these degrees of free- 
dom can be excited easily. Thus, it can be expected that 
they would likewise dominate the transport resonances. 
Another aspect enhancing the effect of dominant CM- 
excitations may arise especially in large quantum dots 
laterally connected to the reservoirs by split gate tunnel 
barriers. This geometry forces the electron to enter at a 
certain point at the edge of the dot. This introduces a 
large dynamical dipole moment in the system built out 
of the incoming electron and the {N — l)-particle state in 
the dot. As known from far-infrared investigations, CM- 
excitations give rise to a strong dipole moment. There- 
fore, one would expect that the overlap of such compound 
state with a CM-excitation will be larger than others, 
thus favoring it for tunneling. 

A basic assumption for the discussion so far has been 
the strict separation of CM- and relative motion which 
is only fulfilled in a truly parabolic confinement poten- 
tial. However, calculations which take into account an- 
harmonicities of the confining potential pl| have shown 
that for realistic conditions ||2^ the coupling between the 
CM-motion and the relative motion is weak, leaving the 
demonstrated mechanism for dominant CM-excitations 
essentially unchanged. 

In summary, the analysis of the overlap matrix ele- 
ments which govern the tunneling rate for an electron 
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traveling through a quantum dot shows that the strong 
correlations present in few-electron dot states are ex- 
tremely important. They strongly reduce the probabil- 
ity of tunneling through channels involving excitations 
of the internal degrees of freedom. This leads to a dom- 
inance of the center-of-mass excitations which are not 
affected by correlation effects. The ccntcr-of-mass exci- 
tation spectrum is identical with the single-particle spec- 
trum of an electron in the parabolic quantum dot, inde- 
pendent of the number of particles. These features give 
an explanation for the experimental observations in the 
non-linear, single-electron tunneling regime, showing rel- 
atively sparse excitation spectra with the characteristic 
single-electron level spacing. 
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FIG. 1. Level scheme of a parabolic quantum dot with 1,2, 
or 3 electrons at S = T. Arrows represent the correspond- 
ing addition energies. Parameters for all figures: Confining 
energy fio = 2 meV, GaAs parameters: dielectric constant 
e = 12.4, effective mass m* = 0.067me. 

FIG. 2. Ground state energy of a parabolic quantum dot 
with 2 and with 3 electrons as a function of the magnetic field 
B. Spin states are indicated: S = Singlet, T = Triplet, D = 
Doublet, Q = Quartet. Arrows mark those field ranges where 
G-G-transitions are forbidden by spin selection rules. 

FIG. 3. Summed values of the overlap matrix elements for 
transitions between the 2-particle ground state in the quan- 
tum dot at -B = and all states of the 3-electron system 
as a function of the transition energy AE. Arrows indicate 
center-of-mass excitations. 
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